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Abstract. We summarize and discuss recent work (Fregeau 2007) that presents the confluence 
of three resuhs suggesting that most Galactic globular clusters are still in the process of core 
contraction, and have not yet reached the thermal equilibrium phase driven by binary scattering 
interactions: that 1) the three clusters that appear to be overabundant in X-ray binaries per 
unit encounter frequency are observationally classified as "core-collapsed," 2) recent numerical 
simulations of cluster evolution with primordial binaries show that structural parameters of clusters 
in the binary-burning phase agree only with "core-collapsed" clusters, and 3) a cluster in the binary- 
burning phase for the last few Gyr should have ~ 5 times more dynamically formed X-ray sources 
than if it were in the core contraction phase for the same time. 
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1. Introduction 

This proceedings article briefly summarizes my contributed talk at the " Population Exp lo- 
sion" meeting, which itself was a brief summary of the work described in iFregeau mm . I 



have taken the opportunity here to inc lude sup p lemen tal material, but flrst refer the reader 
to the detailed, complete discussion in iFregeaul (j2007l ). 



1.1. X-Ray Sources and Cluster Dynamics 

Although X-ray binaries have been known to be o yerabundant in Galact ic globular clusters 
(by unit mass relative to the disk) for over 30 years (jClarklll975l : lKatall975). it is only r ecently 
that clear evidence for their dynamical origin has revealed itself. Pooley et al.l ( 20031 ) found 
that the number of X-ray sources with Lx > 4 x 10^°erg/s, iVx, in a cluster follows a 
clear, nearly linear correlation with the encounter frequency F, a measure of the dynamical 
interaction rate in the cluster. There are a few clear outliers to this trend, however. Terzan 1 
is overabundant in Nx by a factor of ^ 20, NGC 6397 b y a factor of ~ 5, and NGC 7099 by 
a factor of ~ 2 fjCackett et al.l l2006l : iLugger et al.1 120071 ). The common thread among these 
clusters is that they are all classifled observationally as "core-collapsed," while the other 
clusters in the sample are not. A cluster is observationally termed core-collapsed if its surface 
brightness proflle is consistent with a cusp at the limit of observational resolution. Since this 
observational classiflcation is linked with the cluster evolutionary state (as described below), 
it appears that cluster evolution complicates the iVx-F correlation. 



1.2. Understanding Cluster Core Radii 

The evolution of a globular cluster, being a bound self-gravitating system, is very similar 
to that of a star, and comprises three main phases. The cluster first "core contracts" on 
a relaxation timescale. Once the core density becomes large enough for binary scattering 
interactions to begin generating energy, the cluster settles into the "binary-burning" phase. 
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in which the cluster's core properties remain roughly constant with time. Once the binary 
population is exhausted in the core, it collapses via the gravothermal instability, leading to 
extremely high central densities, followed by a series of gravothermal oscillations in which 
the core expands and contracts repeatedly. ( For a graphical r epresentatior i of the three 
main phases of cluster e volution, see Fig ure 1 of iGao et al.l ( 199l[ ). Figure 5 of iFregeau et al. 
or Figure 29.1 of lHeggie fc m^f!2nn3h .) Recent comparison of star cluster evolution 



simulations with observations of cluster structural parameters have shown that the clus- 
ters observationally classified as core-collapsed are in fact most likely in the binary-burning 



phase, while the rest a re most likely still in the core-contraction phase (iHeggie et al.l 12006 
Fregeau fc R^l2007l ). 



1.3. A Refined F 



First introduced by lVerbunt fc HutI f 19871 ). the encounter frequency F is an estimate of the 



current dynamical interaction rate in the cluster, which is assumed to be proportional to the 
current number of observable X-ray sources. Typically, it is approximated as 



C^^int 2 3 / 



where rc is the cluster core radius, pc is the core mass density, and is the core velocity 
dispersion. More accurate approximation s of the interaction rate have been used, including 
numerical integrals over cluster models f Pooley et al.l l2003l ) . but all are estimates of the 



current interaction rate. Recent work has shown that dynamically formed X-ray binaries have 
finite detectable lifetimes (~ 10^ to ~ 10^ yr, depending on binary type), and furthermore 
that there is a lag time o f several Gyr between dyn amical interaction and the binary turning 
on as an X-ray source ( Ivanova et al.l l2006l . l2007l ). Since the binary interaction rate is a 



function of the cluster properties (density, velocity dispersion, etc.), it is clear that the X- 
ray binaries we see in clusters today were formed several Gyr ago, and thus encode the 
recent cluster evolution history in their populations. We thus adopt a refined version of F 
which encodes this history by integrating the interaction rate over the dynamical evolution 
of the cluster. The resulting quantity is the total number of strong dynamical interactions, 
which should be roughly proportional to the currently observable number of X-ray binaries. 
Plugging in the different phases of cluster evolution will yield different estimates for the 
number of sources, and possibly allow one to differentiate among the different phases by 
comparing with the observed number of sources. Using standard assumptions, and adopting 
some results from recent A^-body simulations for the evolution of Tc with time, one finds the 
ratio of the number of interactions for a cluster in the binary-burning phase to one in the 
core contraction phase is 

-^mt,bb _ tx 2.835 , , 

T7 ~ T~r I \ 0.315 7 : \ 0.315 ' 

jVint,cc to ( ^2 1 _ ( *o ] 

\to+9teJ \to+9te+9t:, J 

where t^ is the X-ray source lifetime, t£ is the lag time between strong dynamical interaction 
and X-ray turn-on, and to is the current cluster age. (The gravothermal oscillation phase 
is excluded from the discussion since in this phase the core binary fraction would be much 
smaller than what is observed.) This expression has a minimum of 2.0 and a maximum 
of 17.8 in the range t^ = 10~^-3 Gyr, ti = 1-10 Gyr, for to = 13 Gyr. For the canonical 
values of = 1 Gyr and ti = 3 Gyr with to = 13 Gyr, the value is 5.0. Since the number 
of X-ray sources should scale roughly linearly with the number of interactions, this suggests 
that if a cluster is in the binary-burning phase (and has been for a time tg + t^ to the 
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present), it should have ~ 5 times as many X-ray sources than it would if it were in the core 
contraction phase. The three observationally core-collapsed clusters in the sample previously 
mentioned are overabundant in Nx by factors compatible with the range of values allowed 
by this expression, suggesting that the observationally core-collapsed clusters are indeed in 
the binary-burning phase, while the rest are still in the process of core contraction. 



1.4. Discussion 



Since only ~ 20% of Galactic globular clusters are observationally classified as core-collapsed, 
the conclusion that seems strongly suggested is that most clusters (~ 80%) are currently still 
in the core contraction phase, while those that are core-collapsed are in the binary burning 
phase. This goes counter to the widely held belief that most clusters are currently in the 
binary burning phase, and complicates the many existing studies that have assumed cluster 
core properties that are constant with time, including predictions of blue str aggler popula- 



tions , tidal-capture binaries, and the evolution of the core binary fraction (e.g.. lMapelli et al. 



2004J : iDi Stefano &: Rappaport|[l994 llvanova et al.l [20051 ) . Similarly, this result does away 



with the need for alternate energy sources in cluster cores to explain currently observed 
core radii, including intermediate-mass b lack holes in most clusters, raass lo s s from stellar 
mergers, and ongoing mass segreg ation (jTrenti .2006i : IChatterjee et all l2007l : iMerritt et all 
2004J ). 



2. Meta-Discussion 



A few interesting aspects of this analysis did not fit into the original letter. It has been 
pointed out that an important but neglected factor going into the interaction rate is the 
neutron star retention fraction, since it depends on the cluster escape speed, and since it is 
neutron star -|- stellar binary interactions that predominantly lead to X-ray binaries (Rasio, 
private communication). Since there is no reason to expect that the cluster escape speed 
at the time of neutron star formation does not vary from cluster to cluster, the correlation 
between dynamical interact ion rate and N x should be blurred significantly by this effect. 
We note that the analysis in iFregeaiil ( 20071 ) avoids this issue entirely by comparing a cluster 
only with a version of itself in a different evolutionary phase. 

Another issue is that there is apparently a significant contribution to the X-ray binary 
population from primordial binaries in clusters with small encounter frequencie s. (This is 
likely the reason for the sub-linear dependence of Nx on T in IPooley et al.l (j2003l ) . ) This has 
given rise to the use of the interaction rate per unit cluster mass, 7 = T/Mdus, to help isolate 
the dynamically-formed population in the analysis. The use of 7 has, for e xample, shown 



that a significant frac tion of c a taclys mic variables are dynamically formed (iPooley fc Hut 



20061 ) . The results of iFregeaul (l2007l ) could be further tested in the future with the more 



sensitive 7 as the library of deep X-ray observations of globular clusters grows. 

Finally, we mention that in the near future we plan to test ou r simp le overabundance 
prediction by using the detailed numerical models of llvanova et al.l ( 20051 ) . Although there 
are some hints from that paper that inputting a time-varying core does not noticeably affect 
the final core binary fraction, we expect that the number of visible X-ray binaries, while 
noisy, may be more noticeably affected. 
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